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1.  INTRODUCTION 


The  statistical  description  of  the  geometrical  features  of  the  ocean 
surface  such  as  the  elevation  and  the  slope  is  needed  in  aany  diverse  areas  of 
application  and  has  been  the  goal  of  aany  investigations.  Aside  froa  the 
basic  Interest  of  physical  oceanographers,  the  interest  in  these  features  has 
been  stlaulated  by  the  need  for  reaote  sensing  of  large  areas  of  the  ocean 
froa  satellites  with  various  sensors  including  synthetic  aperture  radar. 
Understanding  of  electromagnetic  scattering  by  the  ocean  surface  is  necessary 
for  this  purpose.  While  theories  abound,  one  of  the  aost  elusive  Ingredients 
in  a  complete  understanding  of  the  problem,  has  been  the  geometry  of  the 
surface,  as  given  by  the  spatial  distribution  of  the  slopes  over  a  sizeable 
area  of  the  ocean  surface. 
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>  Typically,  experimental  measurements  have  been  of  two  kinds:  (a) 
statistical  results  representative  of  a  large  area  and  (b)  a  time  series  of 
measurements  at  a  point.-. The  work  of  Cox  and  Hunk  [1]  and  Stllwell  [2] 
represent  the  first  kind  amd  the  anthology  composed  by  Pierson  and  Stacy  [3] 
is  dominated  by  measurements  of  the  second  kind  where  probes  such  as  wire 
arraysjjr-a-inSer’lilean  measure  the  local  elevation  or  slope  as  a  function  of 
time.  >> Except  for  the  work  reported  by  CUTS7  et  al.  fA)  and  Dobson 
stereo  observations  of  waves,  no  measurements  of  the  Instantaneous  surface 
geometry  over  a  two-dimensional  area  have  been  reported.  Our  own  Interest 
arose  during  a  detailed  analysis  of  the  ability  of  the  Stllwell  method  to 
measure  the  directional  power  spectra  of  ocean  waves. f(| Under  certain 
conditions  of  observation  geometry  and  illumination  conditions,  an  optical 
image  of  the  ocean  surface  can  be  Interpreted  to  obtain  a  good  approximation 
of  the  local  surface  slope  components  along  the  look  azimuth.  If  such  a 
computation  could  be  extended  to  the  slope  component  normal  tooths  look 
direction  and  to  more  general  and  real  illumination  conditions,  a  powerful 
technique  for  measuring  the  ocean  surface  geometry  would  result. 


A  novel  technique  for  reconstructing  the  slopes  of  the  water  surface 
using  three  simultaneous  images  of  the  surface  obtained  with  polarizers 
oriented  at  three  different  angles,  has  been  proposed  by  Krlshnan  [7].  An 
extended  source  of  illumination  such  as  daylight  is  necessary.  The  method 


relics  oa  the  difference  between  the  reflectances  of  the  coaponents  of  light 
polarised  in  the  plane  of  reflection  and  perpendicular  to  the  plane  of 
reflection.  Being  functions  of  the  angle  of  incidence,  the  reflectances  can 
be  used  to  compute  the  plane  and  angle  of  incidence  and  the  slope  of  the 
surface.  It  aust  be  recognised  that  in  the  aeasureaent  of  a  random  process 
such  as  the  ocean  surface  geometry ,  the  only  feasible  check  on  the  results  is 
a  eoaperlson  of  the  statistical  paraaeters  with  the  results  of  other 
aeasureaents  or  through  an  appeal  to  theoretical  arguaents.  However ,  this  new 
method  has  the  unique  characteristic  that  a  detailed  check  of  the  slope 
aeasureaents  of  each  saaple  of  the  random  process  is  feasible.  This  paper 
describes  the  concept  and  the  results  of  its  application  in  a  relatively 
controlled  environment. 

The  basic  theory  related  to  the  concept  and  the  assumptions  and 
limitations  of  the  aethod  are  described  in  Section  2.  Section  3  describes  the 
experlaents  that  were  conducted  at  the  NORDA  test  facility  and  subsequent  film 
processing.  Section  4  describes  the  details  of  the  data  processing  and 
Section  5  shows  exaaple  images  of  the  data  collected  and  the  analysis  results. 


2.  WAVE  SLOPE  RECONSTRUCTION  TECHNIQUE 


2.1  THEORETICAL  FOUNDATION 

The  geometry  specified  in  Figure  2-1,  describes  reflection  at  a  small 

A 

area  of  the  water's  surface  with  the  unit  normal  n  with  zenith  angle  ^  and 
azimuth  v  •  The  angle  of  incidence  (or  reflection)  is  id  ,  and  the  reflected 

A 

ray  r  is  described  by  the  azimuth  a  end  nadir  angle  p  .  The  plane  of 
reflection  at  each  point  on  the  surface  is  defined  by  the  reflected  ray  and 

A 

the  local  surface  normal  n  .  As  shown  in  Figure  2-2,  the  trace  of  the  plane 
of  reflection  in  the  plane  normal  to  the  reflected  ray  (or  the  camera  aperture 
plane)  forms  an  angle  p  with  the  trace  of  the  horizontal  plane,  and  $  is 
determined  by 


»•  ♦  ■  (2.D 

The  angle  of  incidence  u  given  by 

cos  a  -  -sinp  sinp  cos(v-a)  +  cosp  cosp  .  (2.2) 

The  basic  idea  is  to  determine  at  each  location  on  the  surface  the  angle 
of  Incidence  from  the  magnitude  of  the  reflected  intensity  and  the  plane  of 
reflection  from  the  polarization  properties  of  the  light  reflected  from  that 
location.  That  these  in  turn  determine  the  local  surface  normal  and  the 
direction  of  the  incident  ray  at  that  location  can  be  demonstrated  by 
reference  to  Figure  2-3.  The  angle  of  Incidence  defines  a  cone  of  half  angle 

A 

a)  about  the  (known)  direction  r  of  the  reflected  ray  upon  which  surface  normal 

A 

n  must  lie.  Since  the  Incident  and  reflected  rays  and  the  surface  normal  must 

A 

all  lie  in  the  plane  of  reflection,  n  must  be  along  one  of  the  two  lines  of 
intersection  of  this  plane  with  the  cone,  and  the  associated  direction  of  the 
incident  ray  follows  immediately  from  the  laws  of  reflection.  The  more 
probable  direction  for  the  surface  normal  is  assumed  to  be  the  one  lying 
closest  to  the  vertical. 

The  manner  in  which  the  angle  u  and  the  plane  of  reflection  (or  +  )  are 
to  be  determined  will  now  be  demonstrated.  Let  the  sky  radiance  with 
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Figure  2-2  Orientation  of  the  Plane  of  Reflection 
with  Respect  to  the  Horizontal 
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horizontal  polarization  be  denoted  by  (^(A.B),  and  that  with  orthogonal 
polarization  by  N0(A,B),  where  A  la  the  azimuth  and  B  la  the  nadir  angle  of 
the  Incident  ray  (aee  Figure  2-1).  It  can  be  ahown  that  the  plane  of 
reflection  forma  an  angle  <J»  with  the  horizontal  In  the  plane  normal  to  the 
Incident  ray,  where 


coa  4.  -  coa  4  (2.3) 

Thu8,  the  aky  radiance  Incident  on  the  aurface  element  haa  parallel  (p-) 
component 


Np(A,B)  -  Hh(A,B)coaZ«J,  +  No(A,B)sln2<j, 


and  perpendicular  (a-)  component 


Ng(A,B)  -  Nh(A,B)ain2(|,  +  No(A,B)cob2(J,  . 

If  np  and  r)g  are  the  reflectancea  of  water  for  the  p-  and  a-componenta , 
respectively,  the  lrradlance  at  the  image  location  correaponding  to  that 
aurface  element  will  be 


(2.5) 


Ht<«,0)  -  [»pT,p  +  Vsx\a]  (2.6) 

where  F  la  the  f -number  of  the  Imaging  system,  and  subscript  t  denotes  the 
total  lrradlance  (i.e.,  with  no  polarizer). 

If  a  polarizer  oriented  at  angle  y  with  respect  to  the  horizontal  is 
used  In  front  of  the  lens,  the  Image  lrradlance  will  be 

®y(*»0)  "  ^7  {tptVpCO,2<T‘t>  +  N.Vln2(^)] 

+  tsIVp#1"2<T^)  +  *aT,aco,2^T*^^^»  (2.7) 

where,  tp  and  tt  are  the  transmit tances  of  the  polarization  analyzer  parallel 
and  perpendicuar  to  Its  pass  direction,  respectively. 


If  three  simultaneous  Images  of  the  surface  of  water  were  recorded  from 


the  mm  caun  location  with  the  polarizer  pass  directions  oriented  In  the 
horizontal  direction  (h),  the  orthogonal  (o)  direction  and  at  45*  to  the 
horizontal,  the  respective  laage  Irradlances  will  be 

H^a.e)  -  -p-  {tptNphp  cos2*  +  Ngng  »in24> 

+  VVp  *ln2#  +  Vs  cos2*l> 

Ho<a,a>  "  {Wp  8in2*  +  Vs  8i“2*1 

+  VNpnP  co*2*  +  V.8ln^> 

H45(o,B)  -  -Hy  {tptnpnp(l  +  Sln2*)  +  N8n8  (1  -  «ln2*)]  (2.10) 

8F 

+  t8(Npnp(l  -  sln2*>  +  Ngng(l  +  sin2*)] } 
from  which  It  follows  that 

“h  +  Ho  •  ^  (tp  +  '.)(Vp  + 

\  -  Ho  -  ^  (tp  -  e.><Vp '  V.*-2* 

and 

H45  “  {<tp  +  V(Vp  +  Vs>  +  (tp  "  t.)<Vp  "  V.)*ln2*l« 

-  yCH^  +  Ho  +  (^  -  Ho)  tan2*] ,  (2.13) 


(2.11) 

(2.12) 


(2.8) 


(2.9) 


Therefore,  the  knowledge  of  H^,  Hq,  and  H45  determines  tan  2*  and  Equation 

* 

(2.13)  can  be  rearranged  to  give  *  ,  an  estlaate  of  *  ,  by 

A 

♦ 

The  angle  of  Incidence  u  oust  now  be  estlaated  to  complete  the 
solution.  We  use  the  known  functional  relationships  between  w  and 
reflectances  and  np  ,  shown  In  Figure  2-4,  to  coapute  u  • 


■  arctan 


w 


2H 


“5?5h- 


45 


(2.14) 
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The  reflectance  component  Tjg  can  be  found  by  eliminating  N  t|p  from 
Equations  (2.11)  and  (2.12): 


n  - «— 2 -  (H  [t  -t  )cos2^  +  t  +  t  ] 

*  *(t  -t‘)N  cos2*  °  p  8  p  * 

p  8  8 


-tthf(tp-ta)cos2*-tp-tgJ} 


(2.15) 


Note  that  knowledge  of  the  sky  radiance  component  Na  in  the  direction  of  the 
incident  ray  is  necessary  for  solutions  of  this  equation.  However,  the 
direction  of  the  incident  ray  cannot  be  calculated  until  the  surface  normal 
(or  (p,v)],  has  been  calculated,  which  is  the  final  object  of  the 
reconstruction  technique.  An  iterative  solution  with  initial  conditions 
(p*0,  v“0)  nay  be  used  to  overcome  the  above  difficulty.  The  convergence  of 
such  an  iterative  technique  ia  not  always  guaranteed.  He  note  that  for 
uniform  and/or  unpolarlzed  skies  (which  does  not  exist  in  practice)  the  above 
problem  does  not  arise. 

An  alternative  method  of  computing  y  is  as  follows:  A  relationship 
analogous  to  Equation  (2.15)  can  be  derived  as 

2p2 

n  - — 2—2 - (Hh((t  -  t  )cos2*  +  t  +  t  ] 

*<t:-t‘)N,  cos  2*  *  p  8  p  8 

p  a  p 


+  V(tp  -  V0082*  -  cp  -  C.J> 


The  ratio  of  the  two  reflectances  rip  and  is  thus  given  by 


A  ^p  (W'S-V6082*  -  <  W^p+V  N. 

r^  \m  *  <H0-V<tpV  *  *p 


(2.16) 


(2.17) 


For  unpolarlzed  skies. 


hp  •  Ns  or  ! r  "  1  »  < 

and  the  ratio  r  is  independent  of  the  actual  value  of  the  sky  radiance. 


(2.18) 
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Figure  2-4  Reflectance  of  Hater  for  Different  Polarization 
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Once  r  been  computed,  the  known  functional  relationship  between 
h 

id  end  r  (u)  can  be  used  to  find  u  •  although  not  unambiguously*  Figure  2-5 
h 

shows  this  relationship.  We  note  that  each  value  of  r  corresponds  to  two 

values  of  u  ,  one  greater  than  the  Brewster's  angle  at  and  one  less 

c 

than  to  •  Although  this  ambiguity  cannot  be  resolved  completely  by  the 
c 

available  measurements,  it  can  be  eliminated  with  a  great  degree  of  confidence 

by  the  choice  of  appropriate  camera  geometry.  In  the  camera  geometry  employed 

in  our  experiments  (nadir  angle  of  the  optic  axis  -  27*),  u  is  expected  to  be 

less  than  u£(~  53°)  in  almost  all  the  cases.  Therefore,  we  choose  only  that 

value  of  U  which  is  less  than  u  . 

c 

Once  ia  and  b  have  been  determined,  the  surface  slope  can  be  found  as 
follows.  Since  id  and  +  are  the  polar  angle  and  azimuth,  respectively,  of  the 
surface  normal  in  a  spherical  polar  coordinate  system  with  the  reflected  ray 
forming  the  polar  axis,  p  and  v  can  be  determined  by  the  proper 
transformation  of  coordinates  and  are  given  by 

p  -  arccos(slnfa)  sin+  sing  +  coau  cosp)  (2.19) 

v  *  a  +  ATAN2  {-slou  cost,  (sinu  cos$  cosp  -  coau  slop))  .  (2.20) 

This  completes  the  derivation  of  the  surface  slope  from  the  three  polarized 
lrradlances. 

If  the  surface  slope  tilt  and  azimuth  were  known  at  a  given  location,  it 
is  clear  that  the  reflected  lrradiance  can  be  computed  for  a  general  geometry 
and  any  polarizer  pass  direction.  A  comparison  of  the  computed  image  and  an 
actual  image  with  a  polarizer  pass  direction  at  say  -45*  would  then  be  a 
direct  check  on  the  accuracy  of  the  reconstructed  slopes. 

2.2  FEASIBILITY  ADD  LIMITATIONS  OF  THE  TECHNIQUE 

We  now  consider  tha  questions  related  to  the  experimental  feasibllty  and 
the  limitations  of  the  technique  when  applied  to  the  determination  of  the 
slope  of  water  waves.  Some  of  these  are  intrinsic  to  the  procedure,  some  are 

imposed  by  tha  methods  used  for  processing  the  data  and  others  arise  from  the 
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characteristics  of  the  water. 

We  have  already  seen  that  for  uniform  and/or  unpolarized  skies,  the 
analysis  is  more  direct  than  in  the  more  general  case  of  nonuniform  polarized 
skies  where  an  iterative  method  may  be  necessary  for  solution.  While  uniform 
sky  radiance  may  be  uncommon,  unpolarlzed  (but  nonuniform)  sky  radiance  is 
common  especially  when  the  sky  is  hazy  or  overcast.  Even  under  relatively 
clear  conditions,  substantial  parts  of  the  sky  are  unpolarlzed  and  may  be 
useful.  The  ambiguity  in  u  represented  by  Figure  2-5  and  the  observation  that 
the  slopes  of  water  waves  rarely  exceed  0.5  (corresponding  to  a  tilt  angle  p  - 
26.6*)  ,  suggest  the  use  of  a  bias  angle  of  27s,  approximately  one-half  the 
Brewsters'  angle  for  water.  Orienting  the  optic  axis  of  the  imaging  system  at 
a  nadir  angle  of  27*,  limits  the  reflectance  of  water  to  low  values,  reducing 
the  total  amount  of  reflected  light  available  and  the  contrast  in  the 
images.  These  conditions  will  be  particularly  poor  for  the  p-polarlzed 
component  from  those  slopes  which  are  tilted  at  27*  from  the  vertical  away 
from  the  optic  axis.  (Mote  that  the  corresponding  incident  ray  must  be  almost 
horizontal.)  Further,  for  those  slopes  tilted  near  27*  from  the  vertical 
toward  the  optic  ax is,  the  difference  between  the  polarizations  disappears 
because  for  such  slopes  no  unique  plane  of  reflection  exists  as  is  also 
evident  from  the  ambiguity  in  the  value  of  6  f*om  Equation  (2.20)  and 
Figure  2-3.  For  a  given  level  of  noise,  the  errors  in  reconstruction  will 
therefore  be  large  for  slopes  with  tilts  approaching  27*.  However,  with  the 
optic  axis  at  a  nadir  angle  of  27*,  it  is  reasonable  to  expect  that  slopes 
with  tilts  in  the  range  0-15*  can  be  handled.  For  this  bias  angle  and  range 
of  slopes,  obscuration  and  multiple  reflections  do  not  appear  to  be 
significant.  Even  though  the  reflectance  and  the  contrast  will  be  larger  at  a 
bias  angle  of  say  75*,  (on  the  other  side  of  the  minimum  in  Figure  2-5), 
problems  arising  from  obscuration  and  multiple  reflections  argue  against  that 
choice . 

The  fleld-of-vlew  over  which  the  reconstruction  is  feasible  is  limited  by 
several  factors.  Evidently,  for  a  large  fleld-of-vlew,  the  nadir  angle  of  the 
reflected  ray  will  not  be  near  27*  and  the  general  considerations  discussed 

4 

above  will  be  violated.  Zn  addition,  corrections  such  as  cos  z  become 
Important.  Further,  the  transmission  of  optical  systems  are  nonuniform  over 


Figure  2-5  Det ermine t Ion  of  Angle  of  Incidence  from  the 
Ratio  of  the  Reflectances  r  -  n  /n 


large  f islds-of-view .  Finally,  the  behavior  of  polarizers  at  angles  far  from 
normal  Incidence  is  unknown  and  could  only  degrade  the  results. 

So  far  we  had  considered  only  the  reflection  of  the  light  from  the  upper 
hemisphere.  In  reality,  the  light  received  from  the  surface  contains  in 
addition  refracted  upwelling  radiation  and  Equation  (2-6)  should  be  modified 

as 


H(«,p)  -  -V  [»<A,B)t,(«)  +  N'(A\B')  IZlM  ]  (2.21) 

4F  n 

where  n  is  the  refractive  index  of  water  and  N'(A',B')  is  the  upwelling 
radiance  with  incident  direction  specified  by  azimuth  A'  and  polar  angle  B' , 
and  with  the  transmitted  ray  coincident  with  the  refelcted  ray.  The  angles 
(A,B),  (A',B'),  and  u  are  interrelated  by  Snell's  laws.  The  value  of  N  is 
typically  ten  tlsws  the  value  of  H* .  It  is  evident  that  the  effect  of  the 
refracted  light  is  to  raise  the  average  level  of  the  received  light  and  reduce 
the  contrast  in  the  image.  These  effects  nay  be  minimized  by  choosing  to 
operate  in  a  spectral  range  where  upwelling  radiation  is  relatively  small. 
Since  upwelling  radiation  has  a  spectral  peak  in  the  blue-green  region,  a  red 
filter  will  substantially  reduce  the  effects  of  upwelling  radiation.  A  more 
detailed  discussion  can  be  found  in  an  earlier  report  [8]. 

The  success  of  the  technique  depends  on  the  adequacy  of  reflected  light 
within  the  chosen  spectral  region  under  normal  outdoor  illumination  and  the 
slgnal-to-nolse  ratios  achievable  with  the  chosen  detector.  These  questions 
are  examined  fully  in  two  earlier  reports  [7,9].  It  appears  that  high  speed 
infra  red  film  has  adequats  sensitivity  and  that  the  reconstruction  procedure 
is  more  tolerant  of  noise  in  the  value  of  the  slope  tilt  than  the  azimuth. 

Since  the  light  levels  are  already  small,  three  separate  cazwras  must  be 
used  so  that  light  collection  efficiency  is  not  impaired.  This  separation 
necessarily  results  in  errors  due  to  parallax,  l.e.,  the  reflected  rays  (from 
the  same  location  on  the  water  surface)  reaching  the  cameras  are  not 
identical.  The  necessity  of  three  separate  camera  systems  requires  the 
identity  of  their  optical  characteristics,  film  characteristics  and  film 
processing.  Further,  since  pixels  in  the  three  separate  images  corrssponding 
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to  th«  same  surface  location  oust  be  located  for  reconstruction,  there  oust  be 
on  the  surface  at  least  three  fiduciary  narks  so  that  such  a  registration  of 
the  three  inages  can  be  accomplished  (if  internal  registration  is  not 
feasible).  The  presence  of  the  floats  will  cause  some  changes  in  the  surface 
geonetry  and  any  errors  in  the  inage  registration  directly  leads  to  errors  in 
the  slope  values. 

The  moving  surface  of  water  contains  structure  down  to  nillineters 
(capillary  waves)  before  surface  tension  precludes  finer  structure.  The 
spatial  resolution  of  the  slope  estimation  is  therefore  of  interest.  A 
cursory  examination  may  indicate  that  this  would  be  determined  by  the  camera 
and  the  detector  (film  and  microdensitometer).  In  the  present  experiments, 
the  nicrodensl tone ter  aperture  can  be  considered  to  set  the  limit.  However, 
it  should  also  be  noted  that  an  underlying  assumption  is  that  the  reflected 
radiation  averaged  over  thr  (detector)  aperture  is  equal  to  the  radiation 
reflected  by  the  average  slope.  Further,  the  movement  of  the  water  surface 
during  che  exposure  period  is  also  averaged  over  the  (detector)  aperture.  For 
these  reasons,  it  is  necessary  that  the  these  cameras  be  exposed 
simultaneously  and  for  equal  periods  of  exposure,  which  should  be  as  short  as 
possible.  Our  experience  with  the  method  is  not  adequate  to  estimate 
quantitatively  the  errors  introduced  by  these  sources.  Some  of  the  foregoing 
limitations  can  be  removed  by  the  use  of  electronic  recording  of  the  images. 
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3.  THE  EXPERIMENTS  AT  NORDA 


The  experiments  were  conducted  at  the  outdoor  wave  basin  at  the  NORDA 
facility  in  Bay  St.  Louis,  MS  shown  in  Figure  3-1.  Dr.  M-Y  Su  of  NORDA  has 
been  using  the  wave  basin  to  study  non-linear  wave  dynamics  for  -some  time,  by 
mechanically  generating  the  waves  (see  Figure  3-2),  and  following  their 
development  through  their  various  stages  using  capacitance-type  wave  staffs. 
This  instrumentation  Is  severely  limited  in  providing  data  on  the  structure  of 
the  capillary  region  of  the  wave  spectrum  both  by  the  inherent  resolution  and 
by  the  assumptions  and  computations  necessary  to  convert  temporal  spectra  to 
spatial  spectra.  It  would  be  desirable  to  determine  the  latter  by  a  direct 
method  and  remotely,  if  possible. 

Dr.  P.M.  Smith,  also  of  NORDA,  has  been  investigating  microwave  emissions 
from  the  random  wave  surface  in  an  attempt  to  validate  various  models  and  to 
separate  the  contributions  to  the  emlsslvity  of  the  surface  due  to  capillaries 
related  to  wave  instabilities  from  those  directly  induced  by  the  wind.  Drs. 
Daniel  Kwoh  and  Bruce  Lake  of  TRW  have  been  studying  radar  backscatter  from 
the  water  surface  as  a  function  of  wind.  The  interpretation  of  their  results 
will  be  significantly  simplified  if  the  capillary  wave  structure  over 
different  regions  of  the  mechanically  generated  wave-field  can  be 
characterized.  The  optical  measurements  of  wave  slopes  will  provide  data  that 
will  be  useful  to  all  these  researchers  and  will  also  serve  to  establish  a 
unique  new  way  to  provide  sea  surface  truth  for  various  applications. 

Four  identical  Nikon  F2A  SLR  35  mm-cameras  with  Nlkkor  50mm  f/1.4  lenses 
and  MD-3  motor  drives  were  chosen  to  record  the  four  photographs.  The  motor 
drives  could  be  triggered  simultaneously,  from  a  single  remote  switch.  Red 
filters  (Kodak  Wratten  Number  25)  and  polarizers  were  used  in  front  of  the 
lenses.  The  red  filters  were  used  to  block  most  of  the  upwelling  radiation 
dominant  below  600  nm.  The  pass  directions  of  the  polarizers  were  oriented  at 
0,  45,  90  and  -45*  to  the  horizontal.  Kodak  Highspeed  Infrared  film  2481  with 
high  sensitivity  la  the  600-900  as  range  was  used  to  record  the  Images. 

The  cameras  were  attached  to  an  aluminum  frame  with  their  optic  axes  as 
close  as  possible  (to  minimize  parallax  errors).  This  required  the  rotation 


Figure  3-1  (XiCdoor  Wave  Basin  at  NORDA  Facility 


of  two  of  the  cameras  so  that  the  rectangular  Image  area  of  these  two  cameras 
were  oriented  at  90s  to  the  other  two.  The  aluminum  frame  was  attached  to  a 
rotatable  mount  on  an  antenna  tower  which  was  hoisted  erect  after  alignment  r>f 
Che  cameras  (Figures  3-3  and  3-4).  A  television  camera  boreslghted  with  the 
Nikon  cameras  enabled  control  of  the  field-of-vlew  of  the  cameras  from  a 
height  of  70  ft  above  the  water,  with  the  nadir  angle  of  their  optic  axes  at 
about  30°. 

A  plan  view  of  the  experiment  is  shown  in  Figure  3-5.  Four  floats  were 
secured  to  the  bottom  (1  meter  depth)  to  appear  near  the  corners  of  the  FOV. 
Three  of  these  are  needed  for  registering  the  Independent  images  and  the 
fourth  to  check  the  accuracy  of  the  registration.  Four  clearly  resolved  fixed 
objects  whose  positions  are  accurately  known  were  also  in  the  FOV.  These  are 
required  to  calculate  and  check  the  geometry  of  the  optical  systems.  To 
provide  a  calibration  of  the  exposure,  several  gray  cards  and  a  glass  plate 
were  included  in  the  FOV.  To  check  the  simultaneity  of  the  exposure  of  the 
cameras  a  rapidly  rotating  sector  wheel  was  also  used. 

A  large  number  of  experiments  were  done.  Due  to  the  exploratory  nature 
of  these  experiments,  relatively  calm  wind  conditions  and  an  overcast 
unpolarlzed  sky  conditions  were  sought.  A  variety  of  mechanical  waves  were 
photographed.  In  addition,  photographs  were  also  acquired  with  Pan  X  and  Plus 
X  films  with  the  water  heavily  dyed  to  suppress  upwelling  radiation.  A  list 
of  the  available  films  appears  in  Table  3.1. 

Subsequent  to  the  experiments,  the  Hurter-Drif field  characteristic  curves 
for  the  films  used  in  the  experiments  were  generated  by  exposing  strips  of  the 
same  batch  of  films  through  a  calibrated  step  tablet  and  appropriate  filters 
to  unpolarlzed  and  diffused  solar  radiation.  A  schematic  diagram  of  the 
apparatus  used  is  shown  in  Figure  3-6.  A  calibrated  detector  and  an 
oscilloscope  were  used  to  measure  the  energy  density  deposited  on  the  films 
during  these  exposures.  The  four  rolls  of  films  from  each  of  the  experiments 
and  several  calibration  strips  from  the  same  emulsion  were  processed 
simultaneously  by  hand. 

The  processed  films  were  digitized  using  an  Optronics  International  P- 
1000  mlcrodensitometer.  The  H-D  films  were  used  to  convert  photographic 


W.W.  *  wind  waves 
M.W.  =  generated  waves 


.gure  3-6 :  Schematic  Diagram  of  Apparatus  for  Exposing 
Film  Strips  to  Prepare  H-D  Curves 


density  values  to  energy  density*  The  latter  were  used  in  subsequent 
processing . 
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The  calibration  filmstrips  were  also  scanned  and  this  measured  density 
was  related  to  the  energy  density  using  the  data  from  the  detector  output  to 
construct  the  H-D  curves  for  each  film  strip.  The  average  of  several  of  these 
curves  was  used  in  further  data  processing. 


4.  ANALYSIS  PROCEDURE 


Th«  digital  processing  of  the  polarised  vave  luges  can  be  separated  into 
the  following  stages: 

(a)  The  aicrodensl tone ter  provides  an  8-bit  digitisation  of  the 
photograhlc  density  averaged  over  a  square  aperture  on  a  square  grid 
of  sampling  positions*  The  aperture  and  the  sampling  grid  are 
independently  selectable* 

(b)  The  scanned  laages  are  displayed  and  the  coordinates  of  the  floats 
are  deteralned  in  each  of  the  polarised  laages.  These  are  used  to 
calculate  the  translation,  rotation,  and  scaling  needed  to 
superinpose  three  of  the  floats  in  two  of  the  inages  onto  the  third 
(arbitrary)  reference  iaage.  The  registration  accuracy  nay  be 
checked  with  the  reaaining  float.  Using  linear  interpolation,  the 
density  values  of  the  two  superimposed  images  are  calculated  at 
positive  corresponding  to  the  pixel  coordinates  of  the  reference 
image,  thus  generating  a  triplet  of  density  values  corresponding  to 
each  sampled  surface  location. 

(c)  The  data  from  the  Hurter-Dref field  curve  is  now  used  to  convert  the 
photograhlc  density  values  into  energy  density  values.  Using  the 
data  from  the  grey  card  (unpolarlzed  scatterer)  images,  corrections 
are  nude  for  the  exposure  differences.  These  may  also  be  checked 
against  the  images  of  the  glass  plate  (polarizing  dielectric 
reflector). 

(d)  From  measurements  of  the  tower  height,  the  distances  between  several 
fixed  objects  in  the  FOV  and  the  focal  length  of  the  camera,  the 
experimental  geometry  is  reconstructed  so  that  the  correct  azimuth 
and  nadir  angles  (a*{0  of  the  reflected  ray  corresponding  to  each 
pixel  may  be  calculated.  Checks  of  these  computations  agsinst  other 
objects  in  the  FOV  are  also  done  at  this  stage. 

(e)  Using  the  triplet  of  corrected  energy  densities  and  the  reflected  ray 
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geometry,  the  value  of  the  tilt  end  azimuth  of  the  surface  no  real  at 
each  pixel  location  may  be  calculated  as  outlined  in  Section  2. 


(f)  The  reconstructed  surface  normal  results  nay  be  used  to  calculate  the 
polarised  iaage  of  the  fourth  camera  and  several  statistical  results 
on  the  behavior  of  the  two-dimensional  distribution  of  the  slopes. 

The  latter  indicate  histograms  of  the  azimuth  and  tilt  of  surface 
normal,  the  two  components  of  the  slope  along  and  across  the  look 
azimuth  and  the  two  dimensional  spatial  power  spectra  of  the  two 
slope  components. 


5.  OATA  ANALYSIS 


This  section  presents  the  images  and  amplitude  statistics  for  data 
collected  at  Bay  St.  Louis,  Mississipl  on  November  3,  1981  along  with  images 
of  the  sea-slope  constructed  from  the  Images.  Images  from  experiment  IS 
listed  in  Table  3.1  were  used  with  mechanically  generated  waves  and  calm  wind 
conditions.  In  addition  the  difference  between  the  actual  and  constructed 
Irradlance,  which  is  a  check  on  the  accuracy  of  the  reconstructed  slopes  (see 
Section  2)  Is  also  displayed. 

Figures  5-1,  5-2,  5-3,  and  5-4,  Part  a,  are  images  of  the  irradiance  at 
polarizations  of  0*.  45*,  90*,  and  135*  respectively.  Each  figure  has  been 
created  from  photographs,  taken  as  described  in  Section  3,  digitized, 
calibrated,  and  registered.  Part  b  of  each  figure  Is  the  histogram  for  the 
Image  Irr ad lances. 

Figure  5-5  Is  the  image  of  irradiance  at  polarization  angle  135* 
constructed  as  explained  in  Section  2.,  i.e., 

X0°  +  X90°  "  X45° 

Part  a  Is  the  Image.  Part  b  is  the  histogram. 

Figure  5-6  Is  an  Image  of  the  difference  between  the  actual  and  the 
constructed  Irradiance  Image  at  135*  polarization  angle.  Part  a  is  the 
image.  Part  b  Is  the  histogram.  The  amplitude  of  the  error  Is  approximately 
one  fourth  the  Image  amplitude. 

Figure  5-7  Is  an  Image  of  the  surface  slope  angle  u  ,  constructed  from 
the  Images  at  polarization  angles  O’,  45*,  and  90*.  (Figures  5-1,  5-2,  and  5 
3.)  Part  a  Is  the  Image.  Part  b  Is  the  histogram  with  the  angle  measured  In 
radiance.  Figure  5-8  Is  an  Image  of  the  surface  slope  angle  v  ,  constructed 
from  same  images.  Parts  a  and  b  are  the  image  and  Part  c  is  the  histogram. 


Figure  5-la  Irradiance  at  0°  (Horizontal)  Polarization  Angle 


Figure  5-2a  Irradiance  at  45°  Polarization  Angle 


Figure  5-2b  Histogram  of  Irradiance  at  45°  Polarisation  Angle 
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Hlatograa  of  I r rad lance  at  90°  Polarization  Angle 
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Figure  5-6a  Ir rad lance  Difference  Magnitude 
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5-6b  Histograa  of  Ir rad lance  Difference  Magnitude 


Histogram  of  Constructed  Surface  Slope  Angle  y 
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